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Abstract 

Quantum transition and information properties of the Coulomb coupled trion (electron-electron- 
hole) in the double quantum dots under the influence of the time-dependent electric field have been 
studied. Tuning the Hubbard interaction strength amongst the states of electron-hole complexes 
and the parameters of the ac field, two strikingly different kinds of approximate qubit can be con- 
structed within the eight trion states. The similarity and difference between the electron dynamic 
entanglement and overlapping of the wavefunction with the Bell state have been elucidated through 
analyzing these two kinds of qubit. 
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I. INTRODUCTION 



The electronic transport phenomena in the quantum dots have been enriched by the intro- 
duction of the time- dependent driving external fields^, partly due to the sidebands created 
by the harmonic fields, and these sidebands open additional more channels for the tunneling 
of the interacting electrons. More intriguing effect is the dynamic localization produced 
by the harmonic field, which control the electron to locate one of the states permanently 
or in a relatively long duration^^^^^^. The proper choice of the external harmonic 
field parameters can compensate the inadequacy of the artificial engineered nanostructures, 
and provides more possibilities to manipulate the electronic and transport properties of the 
electrons. Since the pioneering works of Loss and Divincenzoii, the interacting electrons in 
the quantum dots have become one of the entities in the realization of solid state quantum 
computation. Besides taking the electron as the quantum information processing candidate, 
exciton&i2 and even charged exciton (trions)^ 1 ^ 1 ^ 1 ^ have been proposed as the candidates 
for the quantum information^^^i. 

The two interacting electrons in coupled quantum dots have been studied by many au- 
thors, where dynamic localization and derealization have been found^^, and possibility 
for the generation of maximum Bell state has been discussed^. Hubbard-like interacting 
excitons in the coupled quantum dots driven by harmonic field have been analyzed and en- 
tanglement has been realized*^. In a single quantum dot, manipulation of the spin of the 
A system consisting of the spin up/down states and trion state is performed by Economou 
and coauthors^ 1 ^, where rotation of spin is realized through an optical way, providing an 
alternative way in controlling spin for quantum information processing. In our work, we 
study the dynamics of the electron-trion X~ in the coupled quantum dots, trion states in 
our model can be formed in either the same quantum dot or spatially-separated dots. The 
manipulation of the transition dynamics among the trion states has been made through 
tuning both the Coulombic interactions among the electron-hole interaction and parameters 
from the harmonic fields. In the Flouquet picture, by resolving the special crossing points in 
quasienergies of the system^iLI^ the resonant tunneling/Rabi flopping within the different 
trion states can be realized. Particular interesting phenomenon lies in the fact that the trion 
is capable of oscillating in the partial states of the all trion states for a long duration of time, 
while the tunneling to the other states becomes forbidden. In these various allowed-forbidden 
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configuration, two different kinds of quantum qubit have been constructed, the overlapping 
of the wavefunction with the maximum Bell state and the concurrence describing the elec- 
tron entanglement have been compared and discussed. The so-called electron entanglement 
dynamic single-slit and double-sli^- have been introduced to describe the similarity and 
difference between the overlapping and entanglement. 

Our work is organized as follows, the model together with the brief description of the 
quasienergies and associated dynamics is presented in section II; The quantum information 
properties including the generation of qubit and dynamic electron entanglement are given 
in the Section III, section IV concludes our work. 

II. THE MODEL AND TUNNELING DYNAMICS 

To describe the electron-trion X~ in the coupled quantum dots, we introduce two quan- 
tum dots system whose energy-level detuning for both the electron and hole is neglected. 
In order for simplicity, The single-level approximation for both the electron and hole is 
adopted, that is to say, the harmonic driving field frequency uo is assumed to be much lower 
than the energy spacing between the single-particle ground state level and the first excited 
state level. If the tunneling rate in the QD systems is much larger than the electron-hole 
recombination rate and the electron-phonon interaction is neglected^ 1 ^, Hubbard on-site 
interaction among electron and hole carriers can be employed, thus the Hamiltonian can be 
written as, 



di t £^, annihilating an electron (hole) i in the i-side (£=\eft (right)) dot. Wi )(J is the hopping 
parameter of electron (hole) between the two dots. U%j^,a an d Vij,a,a denote the intra-dot 
and inter-dot Coulombic interaction among the carriers respectively. nj^ j(J is the number of 
electrons or hole in the left (right) dot. When the QDs are biased by a harmonic field, their 
energy levels will be shifted as follows, ef^ (t) = e° CT =f | cos (iut), where the is the potential 
proportional to the product of amplitude of the harmonic field and spacing between the 




where d\ £ c indicates creating an electron (hole) i in 
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QD. The electron-trion states can be described as follows in terms of the many-particle 
basis d\ t g d\ e a2 dl e \Q) , where i,£, a with the subscripts represent the electron/hole, 
left/right, spin up/down respectively. There are 32 states all together, while only eight 
states survive after taking account of the above assumptions. These eight X~ states can 
be represented as the ket states: \eeh), where the spin- up, spin-down electrons and hole are 
arranged from the left to the right respectively in the ket. By further introducing the number 
and 1 to indicate whether the carriers appear in the left or right QDs, the 8 electron-trion 
states can be written as: |000) , |001) , |010) , |011) , |100) , |101) , |110> , | 111) (denoted as: 
\<f a ),(a = 1,2, ■■■,8)) . Using these basis states, the corresponding Hamiltonian can be 
obtained in a matrix form as, 
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where the gauge transformation exp(— i(2e ea + tho + Uiha)) has been performed to eliminate 
the common term in the diagonal part of the Hamiltonian; the interaction strength \U\\ and 
|Vi| among the carriers in the intra-dot and inter-dot are assumed respectively, without loss 
of generality. In the above effective Hamiltonian, T{t) = (j)cos(ut), and U = —\Ui\ + |Vi|. 
The above system can not be solved in an exact way due to the non-communicativity of 
the Hamiltonian at different times. However Floquet theorem provides a solution clue to 
the time-periodic system : \*& a>m (t)) = exp (— is ajm t) \u a>m (t)). The Floquet state \u a ,rnit)) 
satisfies the following alternative Schrodinger equation, 



u, 



a,m \^a,m/ i 



) 



(3) 



where the quasienergy e Qjm is restricted in the first Brillioun zone [— f , f]- The merit of the 
introduction of the time-dependent external fields lies in that the dynamics can be tuned to 
the regime of the interest. While from the generalized parity symmetry point of view, the 
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subsequent Floquet states can be classified into even parity or odd parity states, the quasi- 
energies belonging to the different and same symmetry may develop into exact crossings 
and avoided crossings respectively^ 1 ^. By choosing the parameters corresponding to these 
special points, the intriguing quantum information properties will be disclosed, which have 
been reported in the section III. 

The changing behavior of the quasienergies and minimum propagator probability p™ m (= 
mm(\c a (tj)\ 2 ),j = 1,2, ■ ■ ■ , N, with = 30T(a = 1,2)) with the driving potential (p has 
been illustrated in Fig.l, where we use the parameters: W e = W 2e = 1.0, Wh = 0.6, u> = 
2.0, the Coulombic interaction difference U has been taken —1. The enhancement of the 
Coulombic interaction U induces asymmetry of quasienergies about the horizontal axis. 
Closer inspection of the Hamiltonian in Eq. ([2]) reveals that the present system can be viewed 
as the two coupled dot arrays with four dots(states) in each array, and hopping term between 
the array is W e (see the Fig. 2). This explains why the quasienegies are split into two parts. 
Also in Fig.l, the survival probabilities P 1 mm and P™ m have been depicted, where the initial 
conditions have been chosen as ci(0) = 1.0 or 02(0) = 1.0 respectively. It is interesting that 
the peaks for p™" and P™ n appear in different exact crossing points of the quasienergies. All 
these show that the states: \ipi) and \ip2) contain different generalized symmetry component 
of Floquet functions^ 1 ^. Another interesting phenomenon is that P™ m have two neighboring 
peaks, forming bistable structures which result from the closely-spacing exact crossings in 
the quasienergies. 

III. THE GENERATION OF QUBIT AND DYNAMIC ELECTRON ENTANGLE- 
MENT 

Rajagopal and Rendell analyzed eight states from triple-qubit, where robust and fragile 
entanglement are classified^. Although there are eights states in our system, these eight 
states from three non-identical particles can not be mapped onto triple-qubit eight states, 
since there exists spatial freedom as well as spin freedom in our system. The states 1 000) 
and | 111) represent that the trion stays in the same QD, while the other states belong to 
the spatially separated trion states. 

Even if there is no obvious qubit in our system, we can tune the system parameters and 
harmonic driving field to create the qubit. Fig. 3 gives transition dynamics among the eight 
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trion states, where parameters used are: U = —20.0, Wh = 0.6, W e = 1.7, u = 2.0, = 24.6. 
The driving potential is chosen to be around the avoided crossing point of the quasienergies. 
From the transition dynamics pattern in the figure, it can be clearly seen that the quantum 
transition nearly takes place among only three states: 1 000) , 1 010) , and 1 100) . Particularly 
interesting phenomenon lies in that the dynamics for |010), and 1 100) is highly equivalent 
in both the amplitude and phase. In this circumstance, 1 000) together with the normalized 
linear combination state: |7/) = -^(|0l0) + 1 100) ) , can be esteemed as a spatial-separated 
qubit^: | 7 ) « ( Cl (t)|000) + c 3 (t)|r/)). 

Inspection of the configuration of the eight trion states indicates that electron dynamic 
entanglement can be realized among the eight trion states from spatial degree of freedom. In 
order to measure the degree of the entanglement, we need to calculate the concurrence^^ 
of the electron-electron density matrix, which is defined as, 



N h 



Q ee := Tr h Q eeh = (lee ®h{Xj\) Qeeh (lee ®\Xj)l 



(4) 



where Tr^ represents the partial trace over the hole basis \Xj)h] ® is the direct product and 
I ee is the 4x4 unit matrix. From the above equation, the electron-electron density matrix 
g ee can be evaluated to be 



| Cl | 2 + |C 2 | 2 C1C3 + C2C4 c x c\ + c 2 c\ cic* 7 + c 2 cl^ 
c 3 c\ + c A c* 2 |c 3 | 2 + |c 4 | 2 c 3 cl + c A c% c 3 c* 7 + c 4 c* 8 



(5) 



C 5 C 1 + C 6 C 2 C 5 C 3 + C 6 C 4 l c 5| 2 + | c 6| 2 C 5 C 7 + C 6 C 8 
1 c 7 c l + C 8 C 2 C 7 C 3 + C 8 C 4 C 7 C 5 + C 8 C 6 l c 7| 2 + | c 8| 2 j 

The other matrix g ee , necessary for the calculation of the concurrence can be constructed 
through the following way, 



ICr y Qee W 
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\ c 7\ 
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(6) 



-(c 3 c* 7 + c 4 c*) c 3 c* 5 + c 4 cl 
y c^ + c^l -(cicg + c 2 c* 6 ) -(cicg + c 2 c\) |ci| 2 + |c 2 | 2 j 

The procedure for the calculation of the concurrence is to find the eigenvalues of the Hermi 



tian matrix «/ y/gg^/g. The concurrence is obtained as C ee = maxjAx — A 2 — A 3 — A 4 , 0}, where 



Aj (i = 1, 2, 3, 4) are the square roots of the eigenvalues of the matrix J tJ~qq^/~q arranged in 
decreasing order— >2£. 

In the case of Fig. 3, let us focus on the electron entanglement of the state 1 010) and 1 100) , 
the non-vanishing Cj(t) are c\(t), c 3 (t), c 5 (t) with the other (i = 2,4,6,7,8) almost 
vanishing, leading to 
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y0 u u j 



whose eigenvalues can be evaluated to be {4|c 3 | 2 |c5| 2 , 0, 0, 0}. Hence, the concurrence 
Cee(t) = 2|c 3 (t)| |c 5 (t)|, which was depicted in the up panel of Fig. 4. The overlapping of 
the state \ip(t)) with the maximum-entangled Bell state \1pBei1) = -^(|010) + 1 100)) , i.e., 
Piif) = \(il>(t)\il>Ben)\ 2 ~ ( I c 3 (t ) + c 5 (£)| 2 )/2 is given in the middle panel of the Fig.4, and 
the bottom panel gives both the curves for comparison. From the figures, it can be clearly 
seen that the occurrence oscillates with time and reach the maximum entanglement when 
the quenching of the trion state 1 000) in the left QD occurs. Particularly, it is remarkable 
that the concurrence is coincident with the overlapping in the present case. This coinci- 
dence can be attributed to the fact that the propagator c 3 (t) is equal to c 5 (t), leading to 
Pi{t) ~ |c3(t) + C5(t)| 2 = 2 1 C3 (t) 1 1 C5 (t) I w C ee {t). The situation here can be pictorially 
described as the electron entanglement dynamic single-slit, where the maximum Bell state 
passes through the qubit (7) ~ (ci(t)|000) + c 3 (t)|?7)). The qubit here only has a dynamic 
single-slit C3(t)\rj) (the other component of qubit ci(£)|000) is blind to the Bell state), which 
is the very reason why overlapping is equivalent to the concurrence. 

The full dynamic entanglement of the electrons taking account of all propagators 
Ci(t), (i = l,2,---,8) is depicted in Fig. 5, where both the concurrence C ee {t) and over- 
lapping pi(t) is given for comparison. The difference between the full concurrence and 
overlapping mainly results from the leaving out the tiny but non- vanishing propagator Ci{t). 
The presence of the tiny propagator cjit) induces the other entanglement: 4j(|000) + |110)). 
The full dynamic concurrence is composed of two sorts of entanglement ^(|010) + 1 100) ) , 
and ^(|000) + |110)), but the former one is much larger than the latter. This explains why 
the full concurrence is slightly larger than the overlapping (see the bottom panel of Fig.5). 
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Does the dynamic double-slit for the electron entanglement exist? As matter of fact, 
besides the above mentioned qubit, there exists other form of qubit which can be consid- 
ered as electron entanglement dynamic double-slit, if we alter both the Coulombic interaction 
strength and harmonic driving field parameters in a proper way. The case is depicted in Fig. 6, 
where the strong driving potential = 40.7 are used, the other parameters used are the same 
as those in Fig.l. The transition dynamics shown in the figure is calculated with the initial 
condition: ( Cl (0), c 2 (0), c 3 (0), c 4 (0), c 5 (0), c 6 (0), c 7 (0), c 8 (0)) = (1/^,0,0,0,0,0,1/^,0). 
It is evident that the transition dynamics for 1 000) , 1 1 10) is virtually the same, making 
transition to 1 001) , 1 1 T 1 ) respectively in a nearly synchronous way. The outcome can be 
attributed to the fact that the strong driving force prevails over the attractive force between 
the electrons and hole, making the hole in one QD tunnel into the other QD. These four 
trion states are able to forming new approximate qubit: \S) ~ (ci(t)\a) + C2(t)\/3)), where 
| a) = ^-(|000) + |110>) and \/3) = ^=(|001> + 1 11 1>). The difference between the qubit \5) 
and the previously defined qubit (7) lies in the different degree of spatial separation of the 
electrons within the two components of qubit. In the previously defined qubit I7), only one 
electron populates in the different QD within the two components of qubit; while all the two 
electrons populate in different QD in the present qubit \S). 

Through taking inner product with the Bell state \/3), the probability P2(t) = \ (j3\^>(t)) \ 2 
for the system to stay in the Bell state \(3) is given in the middle panel of Fig. 7, where 
the up panel gives the concurrence C ee (t) ~ 2|c2(t)||cg(t)|, when we force the propagator 
|ci(i)| 2 and | C7 ) | 2 to vanish . The bottom panel gives both concurrence and overlapping for 
comparison, where the envelope of the overlapping is almost coincident with the concurrence. 
But there is striking difference between the concurrence and overlapping, the pattern in the 
latter figure is irregular, forming many spikes. The reason behind this irregularity comes 
from the tiny phase difference between the two propagators |c2(t)| 2 , and |cg(t)| 2 , which are 
illustrated in the Fig.8. The quantum interference from these two probability terms leads 
to the shaping of the irregular pattern in Fig. 7. Similar to Fig. 7, Fig. 9 gives the overlapping 
P3{t) = \{a\ty(t))\ 2 and concurrence C ee (t) ~ 2|ci(t)||c7(£)| where both C2(t) and c$(t) have 
been set to vanish. The envelope of the overlapping matches the concurrence in a perfect 
way. As matter of fact, both Fig.7 and Fig. 9 correspond to the case of shutting down one 
of the two dynamic slits Ci(t)|a) and C2 (t ) | /3) . It should be pointed out that the state of 
the hole play the role of picking out one of the two the Bell states |a) or \j3). This selecting 
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process corresponds to the overlapping of the wavefunction with state \a) or \(3). 

When both the dynamic slit become transparent, the full concurrence can be approxi- 
mately calculated as 



C ee (t) « Ja+ VBC (8a) 

A = [(| Cl (t)| 2 + |c 2 (t)| 2 ) (|c 7 (t)| 2 + |c 8 (t)| 2 ) + \cl(t)c 7 (t) + c;(t)c 8 (t)| 2 ] (8b) 

B = (\c 7 (t) +c 8 (t)| 2 ) Re(cl(t)c 7 (t) + c*(t)c 8 (t)) (8c) 

C= (| Cl (t) + C2 (t)| 2 )Re(c*(t)c 7 (t)+ C ;(t)c 8 (t)) , (8d) 

where Re represents the real part. The exact numerical value of the full concurrence C ee (t) 
is plotted in the up panel of Fig. 10, the bottom panel of the same figure gives p2(t), P3(t) 
and C ee (t) for comparison, where the full concurrence is coincident with the envelope of ei- 
ther p2{t) or pz(t) in different moment. The concurrence from dynamic double-slit does own 
different evolution behavior than the wavefunction does in quantum mechanics. In contrast 
to the usual double-slit, where the wavefunctions experience the prominent quantum inter- 
ference, the terminology dynamic single-slit and double-slit here are for pictorial description 
since the entanglement interference here bears the striking different behavior than the usual 
double-slit wavefunction interference as illustrated in Fig.5 and Fig. 10. 



IV. CONCLUDING REMARK 



We have investigated the coherent transition dynamics and explored the quantum infor- 
mation properties of the electron-trion X~ confined in two coupled quantum dots driven 
by an ac electric field. The transition dynamics among the eight trion states is controlled 
through the consideration of the symmetry of the Floquet states, which is embodied in the 
avoided crossing and exact crossing in the quasienergies spectrum. The tunneling dynamics 
among the eight trion states also depends sensitively on the difference of Coulombic inter- 
actions between the electrons and hole in the same QD and spatially-separated QD. Tuning 
the motion of the electrons and hole to evolve in partial states of the whole eight trion 
states, the allowed and forbidden trion states within the transition dynamics can be cap- 
tured through which the two different kinds of the qubit can be constructed. From spatial 
degree of freedom, electron entanglement can be produced and the resulting concurrence 
has been given in both analytical and exact numerical results. The similarity and difference 
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between the concurrence and overlapping of the wavefunction with the maximum Bell state 
have been elucidated. Pictorially, one kind of qubit can be conceived as dynamic single- 
slit, while the other kind of qubit can be considered as the the dynamic double-slit for the 
electron entanglement. 
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FIG. 1: Quasienergies e and P™ m {(j)) } (a = 1, 2) plotted as the external ac driving potential 4>, the 
minimum propagator probabilities P™ m ((f)) in a 30 periods time duration experience completely 
different tunneling behavior for a = 1, and 2. 
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FIG. 2: The equivalent schematic picture of the Hamiltonian defined in Eq.([2]), resembling the two 
arrays of single- level dots (states) coupled through the hopping term W e . 
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FIG. 3: The quantum transition dynamics among the eight trion states leading to the formation 
of the first kind of approximate dynamic qubit, the parameters used are declared in the text. 
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FIG. 4: The approximate concurrence calculated from the analytical expression C ee (t) = 
2 1 C3 (t) 1 1 C5 (£) | and exact numerical overlapping of the wavefunction with the maximum Bell state, 
showing the consistency of the two quantity, the parameters used are the same as those in Fig. 3. 
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FIG. 6: The quantum transition dynamics among the eight trion states leading to the formation 
of the second kind of approximate dynamic qubit, the parameters used are declared in the text. 
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FIG. 7: The approximate concurrence calculated from C ee with c\{t) and c?(t) being set to van- 
ish (we weighted Cj(i), (i = 1,3,5,7) by multiplying v /|ci(t)| 2 + | C3 (t) | 2 + | C5 (i) | 2 + | cy(i) | 2 , and 
Ci(t),(i = 2,4,6,8) by multiplying \/\c2(t)\ 2 + l c 4(*)l 2 + l c 6(0l 2 + l c s(*)| 2 )] versus the exact nu- 
merical overlapping of the wavefunction the maximum Bell state, the envelope of the overlapping 
is coincident with the concurrence. 
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FIG. 8: The propagator P2 = \c2(t)\ 2 plotted with blue solid line and Pg = |cs(t)| 2 plotted with red 
dashed line, showing the tiny phase difference between the two propagators probability, leading to 
the spikes of the overlapping in Fig. 7. 
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FIG. 9: The approximate concurrence calculated from C ee with C2(t) and c%(t) being set to van- 
ish (we use the same weight as Fig. 7), and exact numerical overlapping of the wavefunction the 
maximum Bell state, the envelope of the overlapping is coincident with the concurrence. 
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FIG. 10: The full concurrence calculated from the exact numerical result taking account of all 
Ci(t), (i = l,2,---,8) is plotted against the two overlapping, showing the striking difference 
between the full concurrence and individual overlapping in dynamic double-slit for the electron 
entanglement. 
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